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Abstract 

Background: Post-occlusive hyperemia is often used as a paradigm to evaluate vascular reactivity, for example by 
measuring post-ischemic flow-mediated dilation, arterial blood flow or temporally resolved venous blood 
oxygenation (Hb0 2 ). Here we demonstrate the feasibility of a simultaneous measurement of blood flow and Hb0 2 
in the femoral circulation as part of a single procedure. 

Methods: A multi-echo GRE pulse sequence was designed and implemented to collect velocity-encoded 
projections in addition to full-image echoes for field mapping as a means to quantify intravascular magnetic 
susceptibility. The method's feasibility was evaluated at 3T in a small pilot study involving two groups of healthy 
subjects (mean ages 26 ± 1.6 and 59 ± 7.3 years, N = 7 and 5, respectively) in terms of six parameters 
characterizing the time-course of reactive hyperemia and their sensitivity to differentiate age effects. The 
reproducibility was assessed on two of the seven young healthy subjects with three repeated measurements. 

Results: The physiological parameters agree with those obtained with current methods that quantify either 
velocity or Hb0 2 alone. Of the six measures of vascular reactivity, one from each group was significantly different 
in the two subject groups (p < 0.05) even though the study was not powered to detect differences. The mean 
coefficient of variation (CV) from two subjects undergoing repeat scans were approximately 8% for the oximetric 
and the arterial velocimetric parameters in the femoral vein and artery, respectively, considerably below 
intersubject CVs (20 and 35%, for the young and older subject groups, respectively). 

Conclusion: The proposed method is able quantify multiple parameters that may lead to more detailed 
assessment of peripheral vascular reactivity in a single cuff paradigm rather than in separate procedures as required 
previously, thereby improving measurement efficiency and patient comfort. 



Background 

In the United States, cardiovascular diseases (CVD) 
account for more deaths than cancers, chronic obstruc- 
tive pulmonary diseases and accidents combined [1]. 
The burden of CVD in the future is expected to grow 
due to the aging population, prevalence of obesity [2,3] 
and unhealthy lifestyle choices that include smoking [4], 
lack of physical activity [5] and poor nutrition [6,7]; all 
factors predisposing subjects to develop systemic athero- 
sclerosis that can remain silent for decades. The earliest 
stage of pathogenesis of atherosclerosis, the most 
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common cause of CVD, is endothelial dysfunction 
(EDF), which can manifest as impaired reactivity in the 
peripheral vascular bed (e.g. brachial or femoral arteries) 
due to the suppression of nitric oxide (NO) [8] synth- 
esis. In normal vascular function greater demand for 
blood flow in response to a physiological challenge is 
accommodated initially via reduction of microvascular 
resistance. Subsequently, increase in shear stress to the 
vessel wall promotes NO release resulting in vasodila- 
tion (referred to as flow-mediated dilation (FMD)) of 
the conduit artery [9]. Thus, micro- and macrovascular 
function can be assessed during reactive hyperemia. 

Among non-invasive imaging modalities used to quan- 
tify reactive hyperemia in peripheral arteries are near- 
infrared spectroscopy (NIRS) [10-13], single photon 
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emission computed tomography [14], cardiovascular 
magnetic resonance (CMR) [15-18], and ultrasound [19]. 
NIRS and ultrasound (US) are most commonly used 
due to their portability. NIRS is relatively inexpensive, 
has excellent temporal resolution and is less prone to 
artifacts from subject motion. It also allows simulta- 
neous monitoring of relative changes in tissue blood 
flow and oxygenation [12]. However, the method's spa- 
tial resolution is limited and can only provide microvas- 
cular function in relatively superficial tissues. In 
addition, the tissue response to light depends on skin 
color, body fat and muscle layers, leading to large inter- 
subject variations [20]. US measurements of FMD in the 
brachial artery have generated substantial interest in 
clinical research since its introduction [21]. Numerous 
studies have since demonstrated the prognostic value of 
brachial artery FMD [22,23]. However, the methodology 
suffers from substantial inter- and intra-observer varia- 
bility [9,24] in that probe placement or US settings such 
as dynamic range, gain and probe distance can signifi- 
cantly affect diameter measurements [25]. These sources 
of error are of particular concern since the average rela- 
tive change of FMD is on the order of 5% only [26]. 

CMR has also demonstrated its potential for evaluat- 
ing vascular reactivity. CMR flow velocimetry can be 
performed with high spatial and temporal resolution as 
demonstrated in a recent study in which post-occlu- 
sion hyperemia in the femoral artery was evaluated in 
PAD patients and healthy subjects 1.5 T [16]. The 
implementation of the method requires spatially-selec- 
tive rf pulses which are not widely available and migra- 
tion to systems with higher field strength may be 
challenging due to increased field inhomogeneity. 
Alternatively, microvascular function can be assessed 
by quantifying post-occlusive perfusion with arterial 
spin labeling (ASL) [17] in calf muscle. The ASL- 
approach is model-dependent and the temporal resolu- 
tion is limited. 

All of the above methods only target one aspect of the 
vascular reactivity, e.g. hemodynamics, tissue oxygena- 
tion, or FMD, as possible indicators of endothelial dys- 
function. We posit that quantification of several 
physiological parameters may lead to a more sensitive 
assessment for detecting early stage of EDF. As a step 
toward an integrated CMR protocol to detect early 
stages of EDF and to extend a MR procedure developed 
previously based on using oxygen saturation as a 
dynamic tracer [18], we describe an interleaved multi- 
echo gradient-recalled (GRE) sequence that is able of 
temporally resolve venous blood oxygen saturation and 
flow velocity simultaneously during hyperemia. The fea- 
sibility and reproducibility of the proposed method is 
demonstrated in seven young and five older healthy 
subjects. 



Methods 

Integration of MR susceptometry and projection-based 
velocity quantification 

MR susceptometry relies on the quantification of the 
Hb0 2 -dependent susceptibility difference between intra- 
vascular blood and surrounding tissue [27,28]. For the 
remainder of the manuscript Hb0 2 is defined as the 
fractional oxygen saturation expressed in percent. The 
relative susceptibility difference can be derived from a 
phase difference image between two gradient echoes dif- 
fering in echo time. Velocity is computed from a phase 
difference image between interleaved GREs differing in 
velocity-encoding. The basic idea of integrating MR sus- 
ceptometry and velocity quantification is to interleave 
multi-echo GREs with different velocity encoding [27]. 
However, in order to achieve adequate temporal resolu- 
tion to map the time-course of hyperemia, the velocity 
will be quantified with projections [29] while full 2D 
echoes will be acquired for Hb0 2 quantification. 

Pulse Sequence 

The pulse sequence for temporally-resolved measure- 
ment of venous oxygen saturation and flow velocity is 
shown in Figure la. It consists of a fat-suppressed 
multi-echo GRE pulse sequence for collecting velocity- 
encoded projections sampled at TE1 (i.e. in the absence 
of phase-encoding). The same excitation generates two 
additional echoes (TE2 and TE4), which are phase- 
encoded, thereby yielding a full image for field mapping. 
The velocity encoding is toggled after each TR by separ- 
ating the flow-compensating gradient lobes appropri- 
ately (dashed box Figure la); the maximum temporal 
resolution is 2 x TR. Thus, while the pulse sequence 
collects velocity-encoded projections at TE1 the phase- 
encoding gradient is incremented for the remaining 
echoes. In this manner, 64 pairs of velocity-encoded 
projections can be collected while a single 2D field map 
is obtained from TE2 and TE4. The sequence in Figure 
lb collects a 2D reference image at TE1 for removing 
background tissue signal from the velocity-encoded pro- 
jections, as explained in more detail below. In addition, 
data are acquired at TE2 and TE4 for 2D field mapping. 
In this manner, uninterrupted measurement of Hb0 2 is 
ensured. We note that TE2 and TE4 have the same first 
moment in all three directions so that flow-induced 
phase accumulation is minimized for the field maps. 

Subjects 

Seven healthy young (YH, 26 ± 1.6 yrs) and five older 
healthy (OH, 59 ± 7.3 yrs) subjects without prior history 
of cardiovascular disease were recruited to evaluate the 
performance of the interleaved pulse sequence in terms 
of reproducibility and the potential to detect expected 
age-related differences between the two groups with 
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Figure 1 Multi-echo GRE pulse sequence to a) collect velocity-encoded projections (TEl) and full-image echoes, b) Full image is 
collected at TE1 (note phase encodes are repositioned) for removing tissue signal from the projections. The first RF pulse, followed by a spoiler 
gradient, is for lipid suppression. 



parameters of vascular reactivity. To ascertain that all 
subjects were free of peripheral arterial disease the 
ankle-brachial index (ABI) was determined. All subjects 
evaluated had ABI > 0.9. Written informed consent was 
obtained prior to all examinations following an institu- 
tional review board-approved protocol. 

Cuff paradigm 

Ischemia was induced by inflating a blood pressure cuff 
(Aspen Labs A.T.S 1500 Tourniquet System) to 75 
mmHg above the systolic pressure in the upper thigh. 
The location was chosen so as to minimize perturbation 
of the imaging region (approximately 10 cm superior to 
the knee) during cuff inflation or deflation. On ten sub- 
jects (5 YH and 5 OH), the cuff paradigm consisting of 
2 mins baseline, 5 mins occlusion and 6 mins recovery 
was performed to monitor post-ischemic vascular 
reactivity. 

Reproducibility 

Two YH subjects (not taken from the pool of 5 YH) 
participated only in the reproducibility study where an 
abbreviated cuff paradigm (3 mins occlusion and 2 mins 
recovery) was used in order to allow for three repeated 
measurements in a single session to assess reproducibil- 
ity. After each scan the subject was removed from the 
scanner and allowed to rest for 15 mins to ensure full 
recovery before repeating cuff occlusion. 

CMR Protocol 

All imaging was performed on a 3T Siemens Trio. 
Phased-array eight-channel knee coil (Invivo Inc., 
Pewaukee, WI) was used to acquire axial images of 



femoral vessels. From the scout images an appropriate 
readout direction was chosen to avoid overlap between 
the femoral artery and vein along the projection direc- 
tion. For the data acquisition the following imaging 
parameters were used: TE/TR = 5.0/39.1 ms, flip angle 
= 13°, Bandwidth = 520 Hz/pix, voxel size = 1x1x5 
mm 3 , matrix size = 128 x 128. In addition to acquiring 
velocity-encoded projections taken with equal magni- 
tude but opposite-sign first moments, flow- compensated 
projections (not shown in Figure 1) were also acquired 
repeatedly to ensure acquisition of unaliased phase data. 
Thus, three projections offer the flexibility of achieving 
either VENC of 125 or 250 cm/s, as needed, i.e. by com- 
puting the phase difference between acquisitions taken 
with opposite polarity flow-encoding gradients or 
between zero and non-zero first-moment encoded pro- 
jections. The pulse sequence of Figure la was launched 
10 s prior to and was run for 60 s after the cuff release. 
For the remaining recovery period only Hb0 2 was mon- 
itored with the flow-compensated pulse sequence (Fig- 
ure lb). The pulse sequence achieves temporal 
resolutions of 5 s and 117 ms for dynamic oximetry and 
velocity quantification, respectively. 

Data Analysis 
Hb0 2 quantification 

The raw k-space data were processed off-line. A phase 
difference image was computed with the images of 
equal polarity echoes, e.g. TE2 and TE4. The effects of 
low spatial-frequency modulations of static magnetic 
field were minimized by fitting the data, after appropri- 
ate weighting and masking, to a second-order polyno- 
mial [30]. The hemoglobin oxygen saturation, Hb0 2 , in 



Langham and Wehrli Journal of Cardiovascular Magnetic Resonance 201 1, 13:66 
http://www.jcmr-online.eom/content/13/1/66 



Page 4 of 8 



the femoral vein was computed as 

Hb0 2 M 6j_A^|/ATE "I 
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where A(p is the average intravascular phase relative to 
surrounding tissue, ATE = TE4 - TE2, is the inter-echo, 
y is the gyromagnetic ratio of water protons, AXdo - 4tt 
x (0.27 ± O.Q2)ppm [31] is the susceptibility difference 
between fully deoxygenated and fully oxygenated ery- 
throcytes in SI units, 6 is the vessel tilt angle relative to 
the main field B Q and Hct (hematocrit) is the volume 
fraction of erythrocytes in packed blood. 

The time-course of Hb0 2 characterize vascular reac- 
tivity with three parameters: washout time, upslope and 
overshoot, as described previously [18]. The washout 
time refers to the elapsed time to observe oxygen- 
depleted capillary blood at the measurement location 
after cuff deflation (t = 0). Resaturation occurs at a 
mean rate termed upslope and the transient increase in 
blood flow rate leads to an overshoot, an above-baseline 
Hb0 2 level. 

Blood velocity quantification 

In order to isolate the vascular signal from static tissue a 
projection derived from the reference image is sub- 
tracted from the velocity-encoded projections as 
described previously [29]. In brief, the vessels of interest 
are masked out manually from the complex reference 
image and then Fourier transformed to k-space. The 
projection, which contains tissue signal only, is then 
subtracted from the velocity-encoded projections, thus 
isolating blood signal. The subtraction is performed in 
the complex domain to retain the phase information. 
The phase differences between velocity-encoded projec- 
tions are converted to velocity using the relation 
between phase (p and velocity v, as v = (p-VENC/n. The 
spatially averaged velocities of the vessels of interest are 
derived by averaging the velocity along the readout 
direction within the vessel boundaries. Thus, at each 
temporal point, the average velocity represents the velo- 
city averaged over the lumen. 

The following parameters are derived from the time- 
resolved blood flow velocity: time required to reach 
peak average velocity {time-to-peak (TTP)), the duration 
of forward flow (Tff) and maximum change in the flow 
rate relative to baseline (AQ max ), which we define as the 
ratio, v peak> a v g /Vbase, av g - The parameter v base> avg is com- 
puted by averaging the time-resolved velocity over 6-8 
cardiac cycles and v peak> avg is the peak of the temporally 
averaged velocity via 2-second sliding window during 
the course of hyperemia. Computation of AQ max ignores 
the very small change in lumen area that is not detect- 
able at the resolution of 1 mm. This leads to an error 
no greater than 5% based on previous observations 
[16,32] where the lumen diameter was found to increase 
by 3-5% during peak hyperemia. 



Results 

Representative magnitude and phase difference images 
of the femoral artery and vein are shown in Figure 2. 
The phase difference image (Figure 2b) demonstrates 
the susceptibility-induced phase contrast between vein 
and artery. Figure 2c more clearly shows that the rela- 
tive difference in the magnetic susceptibility is negligible 
between tissue and fully oxygenated arterial blood. The 
temporal position of Hb0 2 time-course in femoral vein 
(Figure 2d) is defined with respect to the acquisition 
time of the center of k-space, e.g. 2.5 s, 7.5 s, 12.5 s etc. 
According to error propagation (eqn. 1) the predomi- 
nant contributor to standard error of Hb0 2 is the error 
in Acp; conservatively estimated at ± 3%Hb0 2 . During 
baseline the blood oxygen saturation level is approxi- 
mately 67% (dashed horizontal line). In healthy subjects, 
the observed blood oxygenation at the imaging slice 
typically ranges from 30 to 85% post occlusion, repre- 
senting the saturation level of washed-out oxygen- 
depleted blood and maximum above-baseline value dur- 
ing peak hyperemia, respectively. The range of Hb0 2 
values qualitatively agrees with tissue Hb0 2 measure- 
ments made with near-infrared spectroscopy [12]. Fig- 
ure 3 shows time-resolved velocity data derived from 
projection images of the femoral artery. Figure 3a dis- 
plays time series of projection images, Figure 3b the 
time-course of velocity during hyperemia (expanded in 
Figure 3c). The latter highlights the sharp increase in 
the averaged blood flow velocity and monophasic 
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Figure 2 a) Sample magnitude and b) phase difference image 

for quantifying venous Hb0 2 in the femoral vein of a healthy 

young subject. The region in the dashed box is magnified in panel 

c). d) Representative time-course of Hb0 2 . 
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Figure 3 Post-occlusive blood flow velocity in the femoral artery in a 25 year-old male, a) Magnified view of a velocity image 
reconstructed from velocity-encoded projections during the early phase of hyperemia after the cuff deflation (background cropped). The time 
axis and the readout direction are indicated and each vertical line represents a velocity profile obtained from a pair of velocity-encoded 
projections, b) time-course of blood flow velocity (region of red rectangle magnified in panel c). The duration of forward flow T ff is indicated by 
the double arrow. The magnification in c) shows that each spike corresponds to a systolic peak. 



waveform (absence of retrograde flow) during the early 
phase of hyperemia, presumably due to reduced micro- 
vascular resistance. The velocity waveforms as well as 
the typical peak average velocity agree with previous 
MR quantifications [15,33]. The results are summarized 
in Table 1 where the average value and standard devia- 
tion (SD) for each parameter are reported for each 
group. Uncertainty (determined via error propagation) 



for each parameter is also included under each column 
heading. Based on a two-tailed t-test, assuming unequal 
variances, overshoot and TTP were significantly differ- 
ent (p < 0.05) between YH and OH subjects. We also 
observe that the standard deviations of the parameters 
(all but the overshoot) are greater among the OH, 
which may reflect longer exposure to different types of 
lifestyles. 



Table 1 Summary of physiological parameters derived from dynamic oximetry and velocimetry. 


Age 


Washout 


Upslope 


Overshoot 


AQmax 


Tff 


TTP 


(yrs) 


(± 2.5 sec) 


(± 0.2%HbO 2 /sec) 


(± 3%Hb0 2 ) 


± 0.5 


(± 1 sec) 


(± 1 sec) 








Old (59 ± 7.3 yrs) 








58 


27.5 


0.49 


2 


3.8 


28 


20 


55 


22.5 


0.78 


7 


10.4 


25 


16 


56 


12.5 


1.9 


7 


8.3 


30 


16 


55 


17.5 


0.9 


10 


13.8 


15 


14 


72 


17.5 


1.3 


10 


5.9 


26 


17 


59.2 ± 7.3 


19.5 ± 5.7 


1.07 ± 0.55 


7.2 ± 3.3 


8.44 ± 3.9 


24.8 ± 5.8 


16.6 ± 2.2 








Young (26 ± 1 .6 yrs) 








25 


17.5 


1.2 


13 


9.2 


22 


14 


26 


17.5 


1.8 


25 


12 


25 


14 


26 


12.5 


1.8 


30 


12.5 


30 


11 


25 


12.5 


1.9 


25 


7.7 


20 


11 


29 


12.5 


1.9 


20 


6.65 


25 


14 


26.2 ± 1.6 


14.5 ± 2.7 


1.72 ± 0.29 


22.6 ± 6.4 


9.61 ± 2.6 


24.4 ± 3.8 


12.8 ± 1.64 


*p-value 


0.127 


0.058 


0.003 


0.592 


0.9 


0.017 



^Two-tailed t-test assuming unequal variances. 
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Repeated measurements of Hb0 2 (Figure 4a) and blood 
flow velocity (Figure 4b) are shown for one of the sub- 
jects. The largest variations between scans in the Hb0 2 
time-course typically occur just prior to cuff release due 
to the absence of flow-related enhancement. However, 
this is largely irrelevant since none of the parameters of 
interest are affected by the pre-cuff-release signal course. 
In Figure 4b the blood flow velocity was averaged over 2 
s sliding-window intervals for comparing the velocity 
time-courses of three measurements. The average oxi- 
metric and velocimetric parameters from the three suc- 
cessive scans along with their SD and coefficients of 
variation (CV) are summarized in Table 2. For the repro- 
ducibility study, the average blood flow velocity at rest 



Table 2 Repeat measurements of reactive hyperemia in 
two healthy subjects. 





Washout 
(sec) 


Upslope 
(Hb0 2 /sec) 


Overshoot 
(Hb0 2 ) 


Tff 
(sec) 


TTP 

(sec) 


Vpeak, avg 

(cm/s) 






Healthy 39 yr-old 








Exp 1 


15 


2 


30 


19 


10 


13.5 


Exp 2 


12.5 


1.9 


32 


15 


9 


13.6 


Exp 3 


12.5 


2.2 


30 


17 


9 


12.5 


Avg 


13.3 


2.0 


30.7 


17 


9.3 


13.2 


Std 


1.4 


0.15 


1.15 


2 


0.58 


0.61 


CV (%) 


10.5 


7.5 


3.7 


11.8 


6.2 


4.6 






Healthy 26 yr-old 








Exp 1 


22.5 


1.32 


11 


10 


4.1 


20.6 


Exp 2 


22.5 


1.29 


8 


12 


4.4 


24.6 


Exp 3 


22.5 


1.54 


10 


12 


4 


25.4 


Avg 


22.5 


1.38 


9.7 


11.3 


4.17 


23.5 


Std 


0 


0.136 


1.53 


1.15 


0.21 


2.57 


CV (%) 


0 


9.9 


15.8 


10.2 


5.0 


10.9 



was not quantified (therefore AQ max could not be com- 
puted). Instead, we report peak average velocity, v peak avg , 
(Figure 4b). For both subjects the average CV is approxi- 
mately 8% but the CV was as high as 16% (overshoot in 
26 yr-old subject). However, it should be noted that the 
intersubject SD is about 50% and thus much greater than 
intrasubject SD (Table 1) for all parameters. 

Discussion 

The work shows that simultaneous measurement of 
dynamic oximetric and arterial flow velocity is possible 
for assessing post-occlusive peripheral vascular reactiv- 
ity. In contrast, a sequential measurements would 
require application of the cuff twice, which would not 
be tolerable by a significant number of patients. The 
method returns six parameters, three each related to the 
time-course of venous Hb02, and arterial flow velocity. 
The time-course of Hb0 2 (referred to as dynamic oxi- 
metry) follows the fate of the oxygen-depleted blood in 
the capillary bed upon restoration of flow and its subse- 
quent replacement by normally oxygenated venous 
blood. It thus provides quantitative information on how 
rapidly oxygen- depleted blood is washed out of tissue as 
well as the rate at which tissue is resaturated. In addi- 
tion, the overshoot may reflect the magnitude of hypere- 
mia since tissue has reduced time to extract 0 2 , a time- 
consuming process driven by perfusion, thus greater 
above-baseline Hb0 2 is expected with more pronounced 
hyperemia. The oximetric parameters derived in this 
pilot study are consistent with those of the authors' 
prior work in which a significantly greater overshoot 
was observed in healthy young compared to older sub- 
jects [18]. 
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The velocity profile of the healthy subjects qualita- 
tively agrees with previous "real-time" phase-contrast 
CMR data [15,16] including the duration of forward 
flow (T FF ; 20 - 35 s) and peak average velocity (~ 40 - 
60 cm/s). The unique feature of the "real-time" techni- 
ques [15,16], including the current method, is the ability 
to resolve the velocity waveform for each cardiac cycle 
during hyperemia. In all cases, the blood flow direction 
does not reverse during the early phase of hyperemia 
(Figure 3b) due to significant reduction in microvascular 
resistance. Since vasodilation is generally triggered by 
NO release, TTP and T ff may assess the microvascular 
function's response to an ischemia, i.e. the rate and 
maximum bioavailability of nitric oxide release. 

The projection method has the usual limitation when 
a spatial dimension is sacrificed for improved temporal 
resolution in that velocity is resolved along a single spa- 
tial direction only and thus inherently averaged along 
the projected lumen. Hence, peak velocity cannot be 
quantified but rather an average in the projection direc- 
tion is extracted. Lastly, "misregistration" of femoral ves- 
sels between reference image and projections during 
hyperemia is largely avoided by applying the cuff 10 to 
15 cm from the imaging slice, and since the blood velo- 
city near the vessel wall is nearly zero the error in aver- 
age velocity is expected to be small. 

The proposed method is a component of an ongoing 
development of CMR-based techniques for measuring 
multiple parameters across different vascular territories 
as part of a single examination to study endothelial 
dysfunction. 

Conclusions 

The quantitative imaging method introduced provides a 
series of parameters on peripheral vascular reactivity in 
about fifteen minutes, allowing time for additional quan- 
tification of physiological parameters in a single scan 
session and should be of value to study preclinical cardi- 
ovascular disease. 
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